Strongly correlated states of a small cold atomic cloud from geometric gauge fields 
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Using exact diagonalization for a small system of cold bosonic atoms, we analyze the emergence of 
strongly correlated states in the presence of an artificial magnetic field. This gauge field is generated 
by a laser beam that couples two internal atomic states, and it is related to Berry's geometrical 
phase that emerges when an atom follows adiabatically one of the two eigenstates of the atom-laser 
coupling. Our approach allows us to go beyond the adiabatic approximation, and to characterize 
the generalized Laughlin wave functions that appear in the strong magnetic field limit. 
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I. INTRODUCTION 

Trapped atomic gases provide a unique playground to 
address many-body quantum physics in a very controlled 
way [U [5] . Studies of ultracold atoms submitted to arti- 
ficial gauge fields are particularly interesting in this con- 
text, since they establish a link with the physics of the 
quantum Hall effect. Especially intriguing is the possibil- 
ity of realizing strongly correlated states of the gas, such 
as an atomic analog of the celebrated Laughlin state [3] . 

One way to simulate orbital magnetism is to rotate 
the trapped gas [21 E] ■ One uses in this case the analogy 
between the Coriolis force that appears in the rotating 
frame and the Lorentz force acting on a charged particle 
in a magnetic field. This technique allows one to nucleate 
vortices and observe their ordering in an Abrikosov lat- 
tice [1]. Another promising method takes advantage of 
Berry's geometrical phase that appears when a moving 
atom with multiple internal levels follows adiabatically 
a non-trivial linear combination of these levels [5J. This 
can be achieved in practice by illuminating the gas with 
laser beams that induce a spatially varying coupling be- 
tween atomic internal levels [5] (for a review of recent 
proposals, see e.g. [7]). Recently spectacular experimen- 
tal progress has been made with this technique, leading 
here also to the observation of quantized vortices [S]. 

In this article we focus on the generation of strongly 
correlated states of the atomic gas with geometrical 
gauge fields. We show how Laughlin- type states emerge 
in a small quasi two-dimensional system of trapped 
bosonic atoms when two internal states are coupled by 
a spatially varying laser field. The key point of our ap- 
proach is to go beyond the adiabatic approximation and 
study how the possibility of transitions between the in- 
ternal states modifies the external ground state of the 
gas. We perform exact diagonalization for N — 4 parti- 
cles, in order to analyze the overlap between the exact 
ground state of the system and the Laughlin wave func- 
tion, as a function of the strength of the atom-laser cou- 
pling. We identify a region of parameter space in which 
the ground state, despite having a small overlap with the 



exact Laughlin state, has an interaction energy close to 
zero, a large angular momentum, and a large entropy. 
We show that it can be represented as a Laughlin-like 
state with modified Jastrow factor. 



To better frame our work, let us first emphasize that 
it is well known that the adiabatic accumulation of 
Berry phase induces artificial gauge fields in electroneu- 
tral quantum systems [51 H] . These artificial gauge fields 
of geometric origin are nowadays considered as an impor- 
tant new framework in which strongly correlated quan- 
tum states related to fractional quantum hall physics can 
be engineered in systems of ultracold atoms [7] . The ma- 
jor contribution of the present article is to go beyond 
adiabaticity to show that the main fingerprints of the 
strongly correlated quantum state are preserved in a sig- 
nificantly broad region of the experimental parameter 
space, as demonstrated by detailed numerical calcula- 
tions. This objective is discussed in detail in the frame- 
work of the physics of few-body systems, independently 
of their attainability in the thermodynamic limit, which 
is beyond the scope of the present article. This goal is 
already a realistic one as there are nowadays a number of 
experimental groups able of dealing with small bosonic 
clouds using several techniques [THl E]. These experi- 
mental developments have triggered a number of theo- 
retical proposals focusing on the production of strongly 
correlated quantum states in small atomic clouds [T2HT4] . 



The article is organized in the following way. First, in 
Sec. |Tl] we present the scheme used to generate the arti- 
ficial gauge field together with the formalism employed. 



In Sec. Ill we present our results, discussing in detail 
the properties of the Laughlin-like strongly correlated 
states appearing for different values of the external con- 



trol parameters. In Sec. IV we study the analytical rep- 
resentation of the ground state in the strongly correlated 
Laughlin-like region. Finally, in Sec. |V]we provide some 
conclusions which can be extracted from our work. 
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FIG. 1. (Color online) Schematic illustration of the considered setup, (a) Atoms are trapped in the x-y plane and illuminated 
with a plane wave propagating along the y direction, (b) The energy difference between the two internal states that are coupled 
by the laser field varies linearly along the x direction, (c) Energy eigenvalues of the atom-laser coupling in the rotating wave 
approximation. 



II. THEORETICAL MODEL 

We consider a small quasi two-dimensional ensemble 
of harmonically trapped bosonic atoms in the x-y plane 
interacting with a single laser field treated in a classical 
way. The single-particle Hamiltonian is given by 



P 

Hsp= 2M 



V{r) + H AL , 



(1) 



where M is the atomic mass and V(r) is an external po- 
tential confining the atoms in the plane. Hal includes 
the atom-laser coupling as well as the internal energies. 
In order to minimize the technical aspects of the pro- 
posal, we consider here a very simple laser configuration 
to generate the geometrical gauge field. However, it is 
straightforward to generalize our method to more com- 
plex situations. The laser field is a plane wave with wave 
number k and frequency wl propagating along the di- 
rection y [see Fig. [It a)]. It couples two internal atomic 
states |g) and |e) with a strength that is given by the 
Rabi frequency f?o- 

The atom-laser term comes from the coupling of the 
electric dipole of the atom with the electric field of the 
laser. It can be written as, 

H AL = £ g |g)(g|+i; o |e>(e| 

+ &n o cos(w L t-0)(|e><g| + |g)<e|) , (2) 

with cf> = ky. We neglect the spontaneous emission rate 
of photons from the excited state |e), which is a real- 
istic assumption if the intercombination line of alkali- 
earth or Ytterbium atoms is used. For instance, for 
Ytterbium atoms the state |e) could be taken as the 
first excited level 6 3 Po which has a very long lifetime 
of ~ 10 s [15] . Importantly, we further assume that the 
energies E g = — TiQ,qx/{2w) and E c — Tujja + fi^ox / (2w) 
of the uncoupled internal states vary linearly along x in 
a length scale set by the parameter w, as sketched in 
Fig.[l];b). This can be achieved experimentally either by 
profiting the Zeeman effect, i.e. applying a real magnetic 
field gradient to the system, or by using the a.c. Stark 
shift produced by an extra laser beam with an intensity 
gradient. We suppose that the laser is resonant with the 



atoms in x = 0, i.e. ujl = wa- Using the rotating-wave 
approximation, Hal can be written in the frame rotating 
with the laser frequency and in the {|e), |g)} basis [TS] 
as, 



H AL = 



nn 



e^ sin 9 
— cos 9 



(3) 



where f2 = f2o-\/l + x 2 /w 2 , and tan# = w/x. 

It is convenient to rewrite Eq. ([!]) in the basis of the 
local eigenvectors of Hal, \ipi)r and \ip2)r, associated to 
the eigenvalues hfl/2 and ~Ml/2 : respectively. These 
can be written in the {|e) , \g)} basis as, 



G ( cos 9/2 e 



sin 6»/2 



i<j>/2 



- (-sm.6/2 e 
cos 9/2 e" 



4/2 
-i<p/2 



(4) 



where G = ^ . This particular form of Eq. ([Z ) allows 
us in what follows to obtain a fully symmetric i?22 j see 
Eq. (13). Let us emphasize that the choice of the phase 



factor in front of these eigenstates is nothing but a gauge 
choice for the following. 

The atomic state can be then expressed as 



X {r,t) = ai(r,t) ® |^i) r + a 2 {r,t) ® |V> 2 )r 



(5) 



where a, captures the dynamics of the center of mass and 
\il>i) T of the internal degree of freedom (in the following 
we drop the subindex r in the kets 1^1,2) to simplify the 
notation). Projecting onto the basis •fJV'i}}, and noting 
that 

V r (djipj) = aj(V r ipj) + (V r dj)ijjj , (6) 

the single-particle Hamiltonian, Eq. ([!]), is represented 
by the 2x2 matrix H sp = [H^] acting on the spinor 
[fll(r, t), d2(r, t)]. We find in particular [T71 RlE] 



H 



[p - ejAY 



+ U + V + e. 1 



with ei 



2M ' ' ' ~° 2 ' 

1 and 62 = —1, where we have defined 

A(r) = -ih(4i 2 \V r ip 2 ) 



(7) 



(8) 
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and 



For the chosen gauge they read, 



A{r) = hk 
U(r) = 



y x x 
Aw ' Aw 2^/x r 
\ 2 w 2 



8M(x 2 + w 2 ) 



(10) 
(11) 



We consider atomic clouds extending over distances 
smaller than w. This allows us to expand the matrix 
elements up to second order in x and y. In this 
approximation, we recover the symmetric gauge expres- 
sion A(r) = j^(y,—x) and the artificial magnetic field 
reduces to Bj = Ejhk/(2w) z for an atom in The 
specific choice of the phase factors in Eq. Q is in fact ob- 
tained by imposing as a constraint the symmetric gauge 
at this step. Finally, we fix the external potential V(r) 
such that the total confinement for the spinor component 
a 2 is isotropic with frequency uj±: 

1 -mu J i(x 2 +y 2 ) = U(r) + V(r)- h ^l + ^. (12) 



The Hamiltonian 



Ho 



■ p 2 /2M + p ■ 
, (P + Af 
2M 



A/M 

Muj 2 



Mu]_r 2 /2 



(l-77 2 > 



(13) 



is thus circularly symmetric and its eigenfunctions are the 
Fock-Darwin (FD) functions 4>e jn , with £ and n denoting 
the single-particle angular momentum and the Landau 
level, respectively. The magnetic field strength is charac- 
terized by i] = uj c /2uj±, with uj c = \BA/M = hk/(2Mw) 
the 'cyclotron frequency'. 

The interesting regime for addressing quantum Hall 
physics corresponds to quasi-flat Landau levels, which 
occurs when the magnetic field strength rj is com- 
parable to 1. The energies of the states of the 
Lowest Landau Level (LLL), n = 0, are Eg o = 
huj x [1 + £{l -r]) + {k 2 \l/8) + Xl/{8w 2 )] , where A_l = 

y/h/MLJ_l_. 

Relevant energy scales of the single-particle problem 
are hQ , which characterizes the internal atomic dynam- 
ics, and the recoil energy Er — h 2 k 2 /(2M), which gives 
the scale for the kinetic energy of the atomic center-of- 
mass motion when it absorbs or emits a single photon. 
For hflo 3> £r the adiabatic approximation holds and 
the atoms initially prepared in the internal state j^) 
will remain in this state in the course of their evolution 
[7]. The single-particle Hamiltonian H22, in combination 
with repulsive contact interactions, then leads to quan- 
tum Hall-like physics which has already been extensively 
studied [3] . Our goal here is to consider corrections to the 
adiabatic approximation and to analyze in which respect 
these corrections still allow one to reach strongly corre- 
lated states. This aspect is particularly important from 



an experimental point of view, since the accessible range 
of J7o is limited if one wants to avoid undesired excitation 
of atoms in the sample to higher levels and/or an un- 
wanted laser assisted modification of the atom-atom in- 
teraction. Note that the strength of the atom-laser cou- 
pling, characterized by fio, is distinct from the strength 
of the magnetic field, characterized by rj. Because the 
magnetic field has a geometric origin, rj is independent of 
the atom-laser coupling as long as the adiabatic approx- 
imation is meaningful. 

In the following we consider the situation where hD, Q 
is still relatively large compared to £/r, so that we can 
treat the coupling between the internal subspaces related 
to IV'1,2) in a perturbative manner. In a systematic ex- 
pansion in powers of Qq , the first correction to the adi- 
abatic approximation consists (for the spinor component 
a-i) in replacing H22 by the effective Hamiltonian 16J 



tl 2 2 — tl22 - 



H21H12 



(14) 



The additional term ^21-^12/(^0), which does not com- 
mute with the total angular momentum, is somewhat 
reminiscent of the anisotropic potential that is applied 
to set an atomic cloud in rotation [THl HO] • It is however 
mathematically more involved and physically richer, as 
it includes not only powers of x and y, but also spa- 
tial derivatives with respect to these variables, see Ap- 
pendix [X] for its explicit form. 



III. GROUND STATE PROPERTIES 

The quasi-degeneracy in the LLL can lead to strong 
correlations as the interaction picks a many-body ground 
state for the system. The interaction between the atoms 
is well described by a contact interaction with a coupling 
constant g = \/8ir{a s /l) for the quasi two-dimensional 
confinement. Here a s is the s-wave scattering length and 
I the thickness of the gas in the strongly-confined z di- 
rection. The many-body Hamiltonian then reads 



H 



N 



h: 



eff / 



i=l 



ft 9 \ - 

M 



i<j 



(15) 



Using an algorithm for exact diagonalization within the 
LLL of -H22) we have determined the many-body ground 
state (GS) of the system, providing phase diagrams of 
several relevant average values characterizing the system 
in a broad range of laser couplings, Q , and magnetic 
field strengths, 77. To ensure the validity of the LLL as- 
sumption, we demand that the difference in energy be- 
tween different Landau levels is larger than the kinetic 
energy of any particle in a FD state inside a Landau level. 
In addition, in the full many body problem the interac- 
tion energy per particle is always much smaller than the 
energy difference between adjacent Landau levels. The 
main results are summarized in Figs. [2|[3| a nd |4| and dis- 
cussed in subsections III A III B and III C respectively. 
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In subsection HID we analyze the internal correlations 
in the Laughlin state and in |IIIE we address the impor- 
tant problem of the role of excitations in the Laughlin 
like region. All the calculations are performed for N = 4 
atoms, k = 10/Aj_ 5 and gN = 6. As the perturbation 
H21H12 breaks the rotational symmetry, we cannot carry 
out the exact diagonalization by restricting ourselves to 
a subspace with fixed angular momentum as in standard 
literature. To achieve convergent numerical results, we 
need to include a large number of L subspaces, and this 
number grows as h£l /ER is decreased. In all cases we 
consider all subspaces with < L < L max , where i max 
is chosen to ensure convergency. For dcfiniteness let us 
quote the size of the Hilbert spaces considered in our 
work. For N = 4 we require L max = 28 for most of the 
numerical results reported in the paper, which results in 
a Hilbert space size of 2157. This rapidly growing size of 
the Hilbert space as N is increased, together with our ex- 
plicit interest in providing fine-step phase diagrams vary- 
ing both parameters, 77 and TiQq in a broad region makes 
our full study already computationally very extensive, 
i.e. the calculations presented in this article require on 
the order of 2 weeks on a single 2 GHz processor. 
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FIG. 2. (Color online) Average value of the total angular 
momentum, in units of h, of the ground state for N = 4 
atoms as a function of rj and T1Q.0/ E^. The insets concentrate 
on two different values of hQ,o/En — 40, and 100, respectively. 



A. Angular momentum 

In Fig. [2] we show the expectation value of the total 
angular momentum of the GS as a function of rj and 
Htlo/E-R. For large Slo we recover the step-like struc- 
ture that is well known for rotating bosonic gases, with 
plateaus at L = 0,4,8, and 12 [31 [2_T] with rjoj± play- 
ing the role of the rotating frequency. For an axisym- 
metric potential containing N ^S> 1 bosons it is well 
known that the value L — N corresponds to a sin- 
gle centered vortex, described by the mean-field state 
^ivx = YliLi z i e~'^ z ^ 2Xj - . Here the squared overlap 
between our GS and V^ivx is relatively low (0.47). This 
is due to the small value of N which causes significant 
deviations from the mean-field prediction. The value 
L = N(N — 1) (here L = 12) in the axisymmetric case 
corresponds to the exact Laughlin state, with a filling 
factor 1/2 for any N. For decreasing values of fio the 
transitions between the plateaus become broader and are 
displaced towards smaller values of 77. The Laughlin- 
like region is defined here as the interval of rj fulfilling 
(GS|L|GS) > N(N + 1). 



B. Entropy 

An interesting measure for the correlations in the 
ground state is provided by the one-body entanglement 
entropy [55] defined as, 



S = -Tr \pW In 



(16) 
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FIG. 3. (Color online) Entropy of the ground state for N = 4 
atoms as a function of rj and fifio/^R- The insets concentrate 
on two different values of hQ,o/En — 40, and 100, respectively. 



Here is the one-body density matrix associated to 
the GS wave- function defined as, 



pW(r,r') = (GS|* t (r)*(r')|GS) 



(17) 



where is the field operator, if?(r) — J2e <l>e,oa>e,o, with 
aifi the operator that destroys a particle with angular 
momentum I in the LLL. The natural orbitals, (j>i{r), 
and their corresponding occupations, m, are defined by 
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the eigenvalue problem, 



drp^\r,r')^i{r) =ni<f>i(r'). 



(18) 



The entropy of Eq. ( 16 ) provides information of the 



degree of condensation/fragmentation of the system and 
of the entanglement between one particle and the rest 
of the system. This entropy definition is enough for our 
characterization of the strong correlations of the ground 
state. More detailed studies, such as whether our ground 
state fulfills area laws for the entanglement entropy [23] 
are beyond the scope of the present paper. The entropy 
can be explicitly evaluated as S = — rii lnn^. Thus it 
can be checked that this entropy is zero for a true Bosc- 
Einstein condensate, since all particles occupy the same 
mode (ni = 1,77,; = 0,i > 1). As the system looses con- 
densation, with more than one non-zero eigenvalue, S in- 
creases. For the Laughlin wave-function with N bosons, 
2N — 1 single-particle states are approximately equally 
populated, and the entropy is ~ ln(2iV — 1). The entropy 
S is plotted in Fig.|3]and it presents features that are sim- 
ilar to that of Fig.[2| For a fixed rj, the entropy decreases 
with Qq. For fixed f2 the dependence on r\ exhibits steps, 
similarly to that of (L). The region of (L) = corre- 
sponds to a fairly condensed region with S ~ 0. In the 
one vortex region, corresponding to (L) = N, the con- 
densation is already not complete. This is reflected in 
the abovementioned low value of the squared overlap be- 
tween the GS and ^lvx as well as in the entropy close to 
1. Finally, it gradually increases as we increase 77, and 
reaches its maximum value in the Laughlin-like region, 
77 > 0.93. 



C. Interaction energy 

In Fig. [4] we depict the average interaction energy as a 
function of 77 and Titto/E^. In the inset we also plot for 
L = and L = N = 4 the analytical result expected in 
an axisymmetric potential, E- mt = gN(2N — L — 2)/(8w), 
valid for L = and 2 < L < N [23). The interaction 
energy approaches zero as we increase 77, indicating the 
Laughlin-like nature of the states in the region 77 > 0.93. 

The standard bosonic Laughlin state (at half filling) 
has the analytical form [25 27 



•XZ 1 ,...,Z N )=N'Y[(2 



if 



(19) 



i<3 



where Af is a normalization constant and z = x + iy. It 
is the exact ground state of the system for the contact 
interaction in the adiabatic case|28j. The contribution of 
the interaction to the energy of the system is zero due 
to the zero probability to have two particles at the same 
place. 
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FIG. 4. (Color online) Interaction energy, in units of huj±, 
of the ground state for N = 4 atoms as a function of r\ and 
hQo/E-R. The insets concentrate on two different values of 
HQq/E-r = 40, and 100, respectively. 



D. Internal correlations 

The pair correlation function provides a test for the 
presence of spatial correlations in a system. For the GS 
it is defined as, 



P (2) (r,ro) 



GS 



GS 



(20) 



In Fig. [HJ panels (c,d), p^(r, Fm) is depicted, where f 
is taken as the maximum of the corresponding density, 
depicted in panels (a,b). As seen in the figure, once one 
particle is detected in fo, the other three appear local- 
ized at the remaining three vertices of a rectangle. This 
feature, present also in the exact Laughlin wave function, 
survives both for hfto/En — 40 and 100, even though the 
squared overlap of the ground state with the Laughlin 
differs almost by a factor two, as will be discussed later 
in Sect. |IV| In the adiabatic case, this spatial correlation 
could be inferred from the structure of the analytical ex- 
pression, i.e. the particles tend to avoid each other to 
minimize energy. This is responsible for the particular 
spatial correlation shown in panels (c) and (d). No other 
ground state with L < N(N — 1) exhibits this property. 
A similar phenomenology was found for fermions |29j . 



E. Energy spectrum 

To further characterize the properties of the system 
we discuss the properties of the low energy spectrum, 
and its evolution as we decrease Oo, i.e. increasing the 
non-adiabaticity. In Fig. [6] we show the energy difference 
between the ground state and the first ten excitations as a 
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FIG. 5. (Color online) Density of atoms, panels (a) and (b), 
and pair correlation computed as explained in the text, panels 
(c) and (d), of the ground state for TiQ.q/Er = 40 (a,c) and 
100 (b,d), respectively. The values of rj used are 0.942 and 
0.955 for HQq/Er — 40 and 100, respectively. The solid circle 
marks the position of r*o. The length unit is Ax. 



function of rj. First, let us recall that in the adiabatic case 
the spectrum of the system has already been studied in 
the context of rotating atomic clouds [3], thus our main 
interest here will be to characterize the non-adiabatic 
effects. 

Let us first consider the most symmetric case, Fig. [6] 
(b). In the Laughlin region (rj > 0.952) there are two 
types of lowest excitations: quasi-particle and edge ex- 
citations. For 0.952 < rj < 0.961 the excitation with 
L = N(N — 1) — N marked as A is a quasi-particle type 
state, while for rj > 0.961 the state with L = N(N -1)+1 
marked as Bi is an edge excitation. The tower of edge 
excitations of the system are marked as B n in the fig- 
ure and correspond in the adiabatic case to excitations 
with L = N(N + 1) + n, with n > 0. They are fully 
degenerate in the adiabatic case with a degeneracy given 
by the partition function of n, p(n), defined as the num- 
ber of distinct ways in which n can be written as a sum 
of smaller non-negative integers, i.e. 5 if n = 4 |30j . 
In panel (b) the degeneracy is partly lifted due to the 
slight non-adiabaticity and in panel (a) the condition is 
clearly relaxed. This structure of the edge excitations is 
a fingerprint of the Laughlin state. 

Finally, the maximun energy separation between the 
ground state and its first excitation in the Laughlin- 
region, which in our confined case is both a quasi-particle 
excitation and an edge excitation, increases when de- 
creasing rio ■ It changes from ~ 0.022g huj± for Q,q/Er = 



0.06 
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FIG. 6. (Color online) Energy difference in units of g hui± 
between the first 10 levels of the spectrum and the ground 
state energy as function of rj for hQo/En = 40 (a) and 100 
(b). 
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FIG. 7. (Color online) Squared overlap |(GS|*£)| 2 as a func- 
tion of rj and hQo/En for TV = 4. The dashed line marks 
the region of squared overlap larger than 0.8. The inset de- 
picts the squared overlap for hQ /E-R, — 100 (solid) and 40 
(dashed) as a function of rj. 

100 to - 0.027ghoj ± for n Q /E R = 40. The bulk energy 
difference in the non-adiabatic cases can be estimated by 
linearly extrapolating the segment A to rj = 1, giving 
- 0.18.g hujj_ and - 0.13ghuj ± for Q, /E R = 40 and 100, 
respectively. Thus, by increasing the laser intensity, the 
bulk energy difference approaches the value of the gap 
reported in Ref. [27] for a symmetric and edgeless sys- 
tem. 



IV. ANALYTICAL REPRESENTATION OF THE 
GROUND STATE IN THE LAUGHLIN-LIKE 
REGION 

In this section, we calculate the overlap of the exact 
solutions for the GS in the Laughlin-like region with 
several analytical expressions. To begin with, we cal- 
culate the dependence of the squared overlap | {^J^^c I G^S) | 2 
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of the Laughlin state with the exact GS as a function of 
the magnetic field strength 77 and the atom-laser cou- 
pling ft . The result is plotted in Fig. [7| For large 
Oo (typically > 80 En/h), the adiabatic approximation 
holds (if 22 ~ ^22)'- the overlap between the GS and 
the Laughlin state jumps from a quasi-zero to a large 
(> 0.8) value when the magnetic field strength r\ reaches 
a threshold value. 

For smaller Slo the overlap is much smaller even for 
large r\ (upper left corner of Fig. [7]) . In this case the GS 
must have Jastrow factors that bring the angular momen- 
tum around the value L = N(N — 1) and that suppress 
interactions [Figs. [2] and [2]. Based on these observations, 
we propose an analytical ansatz for this GS of the form 



(21) 



with 



JV 



* £1 = M *£ ' Z i 

i=l 

*£2 = M 2 (*£2 - (*£l|*£ 2 )*£l) 



N 



*£2 = M *£ • ^ , 

i<j 



(22) 



such that we ensure (^cl^a) — and (^ci\^Cj) = 
This ansatz involves components of angular momentum 
L = N(N - 1) and L = N(N - 1) + 2, and zero in- 
teraction energy. The coefficients a, /3 and 7 are given 
by the projections of the exact GS onto $£, $£i and 
$£ 2 respectively. In Fig. [8] (a,b) we present the squared 
overlaps ^Laughlin and Pl ± between the exact GS wave 
function and the functions 5 , £ and ^o, respectively. We 
restrict our study to the Laughlin-region. We also plot 
the weights of the angular momentum subspaces in the 
GS, Pl=n(n-i) and Pl=n(n-x)+2- 

The first result of our numerical analysis is that Pc 2 
is negligible (< 0.005) over the whole range of Fig. [8] 
Then, we note that the relations Pl=n(n-i) ~ Pc and 
Pl=n(n-i)+2 ~ Pci hold over this range. This implies 
that the deviation with respect to the adiabatic approxi- 
mation mostly increases the weight of the *&ci component 
in the GS. For small values of MIq/E^, the squared over- 
lap with the proposed ansatz reaches values of ~ 0.85, 
with the weight of <i/£ and ^>ci being of comparable size. 
As MIq/Eh increases above 80, the GS is very well rep- 
resented by Eq. (fl9|), as already explained. Considering 



different particle numbers from N = 3 to N = 5, we al- 
ways find a very similar behavior and thus conclude that 
Eq. (21) quite generally provides a good representation 
of the GS in the Laughlin-like region^] 



1 Note that the Laughlin-like region decreases notably in size as 
N is increased. 
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FIG. 8. (Color online) (upper panels) Squared overlaps be- 
tween the GS of the system and the Laughlin wave func- 
tion, PLaug = |(GS|*£)| 2 (solid squares) and ^ci, Ph\ = 
|(GS|*£i)| 2 (triangles). The sum of both is depicted as solid 
diamonds. The solid and dashed lines correspond to the 
weights Pl=n(n-i) and Pl=n(n-i)+2 in the GS, respectively. 
Panel (a) corresponds to hQo/En = 40 and (b) to 100. The 
shaded band marks the region of squared overlap larger than 
0.95. (lower panels) Energy difference in units of hoj± be- 
tween the first 10 levels of the spectrum and the ground state 
energy as function of r\ for HQq/E-r, = 40 (c) and 100 (d). 



SUMMARY AND CONCLUSIONS 



In conclusion we have performed exact diagonalization 
to analyze the ground state of a small cloud of bosonic 
atoms subjected to an artificial gauge field. Our ap- 
proach allowed us to explore both the regime of very 
large atom-laser coupling, where the adiabatic approx- 
imation is valid, and the case of intermediate coupling 
strengths. In the first case we recovered the known re- 
sults for a single component gas in an axisymmetric po- 
tential. The second case is crucial for practical imple- 
mentations because it requires less light intensity on the 
atoms, which decreases the residual heating due to pho- 
ton scattering. In this case we have identified a regime 
where a strongly correlated ground state emerges, which 
shares many similarities with the Laughlin state in terms 
of angular momentum, energy, internal spatial correla- 
tions and lowest excitations, although the overlap be- 
tween the two remains small. Importantly, a reduction 
of the laser intensity shifts the region where Laughlin-like 
states exist to lower values of the effective magnetic field, 
thus departing from the instability region, rj > 1. We 
have also proposed an ansatz that represents the ground 
state quite accurately for a region of the parameter space. 
Finally, let us emphasize that the properties analyzed in 
this article are measurable quantities, as is the case of 
the expected value of the angular momentum, the pair 
correlation distribution and excitation spectrum. 
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Appendix A: Explicit form of the H21H12 term 

We provide here the explicit expression for the term 
H21H12 appearing in the perturbatively derived Hamil- 



tonian H^- As explained in the text we consider up to 
quadratic terms in x and y. The explicit expression then 
reads 



H21H12 



ft 4 



2x 2 h 4 



k 2 x 2 h 4 



\4M 2 w 4 
ikxyh 4 



M 2 w 6 16M 2 w 4 
k 2 y 2 h 4 



4M 2 w 5 64M 2 w 4 



ikxh 



ik 3 xh 4 



AM 2 w 3 8M 2 w 



xTf 
M 2 w 4 

k 2 h 4 



AM 2 
h 4 
AM 2 i 



+ 



ikyh 
8MV 

k 2 x 2 h 4 

AM 2 w 2 J 
x 2 h 4 
2M 2 w 4 



0: 



k 4 x 2 h 4 

MM 2 w 2 
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